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Abstract 



We carry out a detailed analysis of fluxes of cosmic ray antiprotons, positrons, electrons 
and photons to be expected from the annihilation of relic heavy neutrinos in the galactic 
halo. The spectra of particles are evaluated by Monte Carlo simulation. The results of 
calculations show that the investigation of cosmic ray positron spectra at high energies 
could provide a distinctive signal for annihilation of very heavy neutrinos in the Galaxy 
and give an important information on parameters of dark matter particles. 



1 Introduction 



The agreement that more or less exists between the observed and predicted abundances 
of light elements such as D, 3 He, 4 He and 1 Li provides a confirmation of the standard 
cosmology and leads to the determination of the baryon density in the range |1[ 

p B = (1.7 - 4.1) ■ 1(T 31 g/cm 3 . (1) 



Because the critical density 
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p c = 3H 2 /8nG = 1.88 h 2 ■ 10~ 29 g/cm 3 



(2) 



depends on the Hubble constant H (h is the Hubble constant in units of 100 km s 1 Mpc : ) 
the present ratio of the baryon density of the Universe to its critical density also depends on 

tt B = Pb / Pc = (2.2 - 9.0) 10- 2 h' 2 . (3) 

Therefore for a wide range of the Hubble constant values, h ~ (0.4 — 1), baryons contribute 
(1 — 15)% of the critical density of the Universe. Since the contribution of luminous matter 

Q lm « S-HT 3 /*- 1 (4) 

is much less than 1% of the closure density, most baryons in the Universe should be dark 
either in the form of hot gas in clusters of galaxies or in the form of dark stars. 

There are increasing observational evidences on a variety of scales favouring an average 
density of matter in the Universe which is significantly greater than 15% of the critical 
density possibly contributed by baryons. The non-Keplerian character of rotation curves 
around spiral galaxies indicates that they are at least ten times more massive than the sum 
of the stars which are being seen JJ. The velocity dispersion of stars and X-ray emission 
in elliptical galaxies show that these objects are also dominated by a dark component. At 
the scale of clusters of galaxies, independent sets of measurements point at an even greater 
amount of the dark matter. If this is the case, most of the matter in the Universe must 
be nonbaryonic. The evidence for dark matter in galaxies is collected in Ref.|3[], and the 
evidence for dark matter in the Universe is summarized in Ref.0]. 

Nonbaryonic dark matter can be divided into two broad classes on the basis of its clustering 
properties — the hot dark matter being relativistic at galaxy formation time and the cold 
one being slowly moving at that epoch. The most favoured hot dark matter candidate is 
a light neutrino with a mass ~ 10 eV. Since light neutrinos are relativistic until relatively 
recent epochs, they would free-stream out off and damp out density pertubations up to the 
scale of galaxy clusters. However, the results from COBE show that the observed large scale 
structure cannot be explained with purely hot dark matter ||. 

The cold dark matter on the another hand clusters at all scales, since it must have negligible 
free-streaming length. In the cold dark matter models, the structure generally forms hierar- 
chically, with smaller clumps merging to form the larger ones. For this reason, theorists used 
to prefer the cold dark matter scenario over the hot one. But it is worth to note the recently 
surging popularity of the mixed scenario: a combination of 70% cold dark matter and 30% 
hot one produces a favourable spectrum of density pertubations and leads to consistency 
between the inflationary predictions and some dynamical estimations ||. Theoretically, the 
favourite candidates for cold dark matter are the axion, an ultra-light pseudoscalar particle, 
and weakly interacting massive particles (WIMPS), with masses generally of the order of 
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100 GeV |7], §]. The most attractive WIMPS candidate is the neutralino — a supesymmetric 
fermionic partner of the Standard Model bosons. For a review of the supersymmetric dark 
matter see ||. In some models heavy neutrinos play also an important role in cosmology 
contributing part of the closure density |TI], |TTj, |T2] . 

The nature of the dark matter is certainly one of the most important questions facing both 
particle physics and cosmology. An indirect way to detect such dark matter particles is to 
study their signatures in cosmic rays. 



In the present paper we carry out a detailed analysis of the influence of the annihilation of 
very heavy neutrinos (m v > 45 GeV) on cosmic ray production in the Galaxy. There exist 



a number of scenarios O, ITJ in which such very heavy neutrinos are expected to occur. To 



be explicit, we consider the standard electroweak model, including, however, one additiona 
fermion family. The heavy neutrino v and the heavy charged lepton L form a standard 
SU(2)l doublet. In order to ensure the stability of the heavy neutrino v we assume that 
m u < m L and that the heavy neutrino is a Dirac neutrino. 



The organization of the paper is the following. In Sec. 2 we discuss the residual relic number 
density of heavy neutrinos in the Universe and in the Galaxy. In Sec. 3 we present various 
detailes of our calculations. In Sec. 4 we discuss the results for the spectra of electrons, 
positrons, antiprotons and photons resulting from neutrino annihilation in the galactic halo 
and discuss signatures for the detection of heavy neutrino in cosmic rays. In Sec. 5 we 
summarize the results of our analysis of cosmic ray spectra from neutrino annihilation. In 
the Appendix we present formulae for matrix elements of the reactions considered in the 
paper. 



2 Relic number density of heavy neutrinos 



It is known that the available experimental results do not preclude the existence of heavy 
Dirac neutrinos with masses m v > rriz/2 |15| . In the early Universe at high temperatures 
(T » m u ) such heavy neutrinos (if they exist) should be in thermal equilibrium with other 
species of particles. As the temperature in the Universe drops, heavy neutrinos become 
nonrelativistic at T ~ m v and their equilibrium concentration is proportional to n u ~ 
exp(—m u /T) |frB[. However, in the further expansion of the Universe, as the temperature 
drops below the freeze-out value T/, the weak interaction processes become too slow to keep 
neutrinos in equilibrium with other particles. As a consequence, the number density of heavy 
neutrinos fails to follow the equilibrium concentration (the exponential drop of the number 
density becomes much slower) reaching at present [PL8| : 



n v {T) 



2 • 10" 18 M n 



m v VreloMl 



40 + In 



g s m u r2 

— Jj-VreiaM 



n 7 (T), 



(5) 



where n 7 (T) = 0.24T 3 is the equilibrium photon number density, T = 2.7 K is the present 
photon temperature, M p is the proton mass, v re i is the relative vv velocity in units of 
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the velocity of light. The bar over quantity v re i<j means the averaging over the velocity 
distribution of the heavy neutrino. 



The freeze-out temperature Tf can be found from the equation [16]: 



In 



m,. 



T 



f 



m,. 



-— exp 



T 
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40 + In 



g s m v 2 

— — V rd aM 

g*M„ 



(6) 



In the above expressions g s is the number of particle spin states (for photons and massive 
fermions g s = 2). g*(T) = Nt, os + ^Nf erm is the number of effective degrees of freedom 
at given temperature (see JTB]]), Nf, os is the number of bosons and Nf erm is the number of 
fermions (without heavy neutrinos) being in the thermal equilibrium at given temperature. 

It was assumed in Eq.(^jj) that neutrinos are nonrelativistic at freeze-out, Tf « m u . The 
numerical solution for Eq.(|6]) gives Tf = 3.90 GeV and Tf = 10.7 GeV for the heavy 
neutrino masses equal to m„ = 100 GeV and m v = 300 GeV, respectively. Thus, the 
non-relativistic approximation is valid. Therefore we can neglect the dependence of the 
annihilation cross section on velocity of the neutrino and in further calculations we take 
c(y / s) = a{2m v ). In this case we can replace v re i a « n^a(2 m v ). Because v1„, = 2v 2 for 



definitness in further calculations we can put v re i = \/2v, where v 
average velocity of heavy neutrinos in Galactic halo. 



Because v^ el 

300 km ■ s^ 1 is the 



We have to note that Eq.flg) has been obtained using the analytical solution 0] for the 
expansion equation of hot equilibrium plasma. Such an approach gives correct order of 
magnitude, but, in principle, the equation has to be solved numerically to obtain more 
precise results. In the present paper we do not do this because our main aim is to show 
influence of the hadronization of final state quarks on the e ± , pp and 7 fluxes as well as to 
correctly take into account energy losses for charged final state particles travelled through 
the Galaxy. 



In the general case the following processes could lead to the annihilation of heavy neutrinos 
in the early Universe: 



v v 



ff 

w+ w- 

H°H° 
Z° Z° 
Z°H° 



s channel Z°, H° exchange 

s channel Z°, H° and t channel heavy lepton L exchanges 
s channel H° and t channel heavy neutrino v exchanges 
s channel H° and t channel heavy neutrino v exchanges 
s channel Z° and t channel heavy neutrino v exchanges 



The annihilation of heavy neutrinos into a final state containing Z° boson or/and H° meson 
are suppressed in comparison with W + W~ final state and will be neglected in the following 
consideration. 

Qualitatively, the dependence of the present neutrino number density on its mass m v can be 
easily understood if one takes into account that neutrinos at freeze-out are nonrelativistic 
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and a depends only on the mass of the heavy neutrino. The rate of "burning-out" of heavy 
neutrinos is proportional to the rate of collisions N ~ van 2 u . Using Eq.(|5|) we obtain 
N ~ l/o- at T ~ Tf. 

Thus the concentration of the heavy neutrinos at present, roughly speaking, depends in- 
versely on its annihilation cross section. Therefore, in the region m v ~ the neutrino 
number density is extremely small as a result of the huge value of the annihilation cross sec- 
tion at the Z° boson pole. When the neutrino mass increases the cross section for neutrino 
annihilation into fermions drops and this leads to an increase of the neutrino number den- 
sity, which reaches about 3 • 10~ 10 cm" 3 at m v ~ 100 GeV. At m v > my/ the annihilation 
channel into W + W~ opens and gradually becomes the dominant one since its cross section 
grows like ml |T7| and the present v number density drops again. As it was shown in 



the neutrino number density at m u > 2 m\y is given, approximately, by 

2 • 10" 9 



n,, ~ cm 3 . (7) 



(m u /m w ) \{m v /m w ) +4 

Above we discussed the neutrino number density in the Universe. However, at the stage of 
the formation of the Galaxy the motion of heavy neutrinos in the nonstatic gravitational 
field of ordinary matter, which contracts as a results of energy dissipation via radiation, 
provides an effective mechanism of energy dissipation for neutrinos |19| . As a consequence, 
the contracting ordinary matter induces the collapse of the neutrino gas and leads to the 
following significant increase in the neutrino density in the Galaxy |TJ, 19]: 




where, p G ~ 5 • 10 24 g/cm 3 is the average density of matter in the Galaxy and p 
4 • 10~ 31 g/cm~ 3 is the density of matter in the Universe. 



3 Fluxes of particles from the heavy neutrino annihi- 
lation in the halo 



The annihilation of dark matter particles in the galactic halo could produce detectable 
fluxes of elementary particles such as electrons, positrons, antiprotons and photons. For 
the case of supersymmetric dark matter particles a number of calculations of cosmic ray 
fluxes have been made (see, e.g., |2(], [H], ^2|). In this Section we present the results of 
our numerical calculations for fluxes of particles which appear either directly in the heavy 
neutrino annihilation or as secondary products of the fragmentation of quarks and gluons 
and from the decays of unstable particles. 



Let us consider first the fluxes of electrons and positrons from heavy neutrino annihilation. 
In this case the differential spectrum of electrons (positrons) observable at the Earth is 



defined by [[19 
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where r e is the confinement time for electrons (positrons) in the galactic halo equal to 
r e pa 10 7 years and v is the interstellar e~(e + ) velocity. The quantity dn e /dE is the energy 
spectrum of electrons produced in the annihilation of the heavy neutrino and antineutrino 
and normalized on one collision; n u is the present density of relic heavy neutrinos in the 
galactic halo. 

In our calculations we assume that r e pa 10 7 years, but we have to note that the uncertainties 
here are substantial and can reach an order of magnitude. Moreover, r e can, in principle, 
depend on the initial energy of electron. Additional complications arise since the electrons 
(positrons) trapped in the galactic halo lose energy due to various processes. The energy 
losses for electrons and positrons propagating through the Galaxy halo matter are being 
defined by pi M: 



dE r i 

— = ~ [2.7 n H + 7.3n H E + 1.02 {U rad + U mag ) E 2 ] x 1(T 16 GeV s~ l 



(10) 



where U ra d is the ambient photon energy density in eV cm~ 3 and E is the e~(e + ) energy 
in GeV. The four terms in ( [T0|) represent the energy losses due to ionization of the H 
atoms, bremsstrahlung, inverse Compton scattering and synchrotron radiation in the galac- 
tic magnetic field, respectively. The electrons (positrons) spend most of their time, and 
lose most of their energy in the halo outside of the disc, where tin ~ 0.01 atoms cm~ 3 , 
Urad ~ 0.25 eV cm~ 3 and U mag ~ 0.05 eV cm~ 3 , so that inverse Compton scattering is 
dominant for E > 1 GeV. The e~ and e + spectra produced in neutrino annihilation must 
be corrected for these energy losses before the comparison with observations. 

To calculate the particle fluxes from annihilation of very heavy neutrinos in the Galaxy we 
used the Monte Carlo approach. Such an approach allows to take into account the production 
of the particles of interest (e + , e~ , 7, p, p) by hadronization of initial quarks or decays of 
heavier states. In this approach it is also easy to take into account the energy losses of final 
state electrons and positrons in the Galaxy. 

Because we consider large masses of the heavy neutrino the resulting energies of final particles 
are also large, E > 1 GeV. In this case the dominant contribution comes from the two 
last terms in Eg-flHTf), which are proportional to E 2 , so that, dE = —aE 2 dt, where a = 
1.02 • 10~ 16 {U ra d + U mag ). Therefore, to simulate the broadening effect of energy losses we 
can use the distribution followed from the simple, spatially homogeneous model 



which leads to: 



m*) , (id 



E = E ° . (12) 



1 - 1.02 (U r ad + U mag ) 10- 16 £ r e ln£' 



6 



where E is the final energy of particle, E is its initial energy and £ is the random variable, 
with uniform distribution in the (0, 1). 



As it was mentioned in Sec.0 we restricted our consideration to the following two processes: 

u + u ^ f + f, , s 

v + v -> W + + W~, { ' 

where / denotes a lepton or a quark. The expressions for the differential cross section and 
matrix elements are given explicitly in the Appendix. 

The simulation of the annihilation processes were perfomed with the package PYTHIA 5.7 



27] with suitable modifications to include the fourth generation of fermions. To perform 



the hadronization and decays of final states we used the package JETSET 7.4 P7| . This 
package is well tuned to the experimental data and gives a satisfactory description of final 
state multiplicities and particle distributions. We assumed that there is no mixing between 
the three light generations and the fourth one. For the hadroniztion of light quarks and 
gluons we used the symmetric LUND fragmentation functions and for heavy quarks the 
Peterson function. We also took into account the final state radiation in the annihilation 
reaction. This leads to a drastical increase of the spectra in the very low energy region. We 
have taken into account the decays of w, n and fi, which particles are usually considered to be 
stable. To avoid unnecessary complications we chose very high masses for the heavy lepton 
and Higgs = 1000 GeV, m H = 10000 GeV, because, as we have already mentioned 
previously, the Higgs contribution is important only in the resonance region m v ph tuh/2. 
We performed the simulation for two mass values of the heavy neutrino: m v = 100 GeV and 
m v = 300 GeV assuming for the mass of the t quark m t = 170 GeV. The other parameters 
have been left equal to their default values used in the PYTHIA 5.7/ JETSET 7.4 package 
(see for detailes 



The results of our calculations for particle spectra are presented in Figs||, |] and for total 
cross sections in Table. 1. 



4 Discussion of numerical results 



The results of the numerical simulation of the electron flux are presented in Fig.|3](a) and in 
Fig.f|(a) for m v = 100 GeV and m v = 300 GeV, respectively, together with an approxima- 



tion of the observable electron spectrum |23[ of the type 



J™P = 0.07^ 3 - 3 crrT 2 sr^s^GeV- 1 , 
where E is e~(e + ) energy in GeV. 



(14) 



Without energy losses |T8[ we could expect a very narrow peak at E e = m u /2 with width 
mv/c well above the background from the annihilation process v v —>■ e + e~ (since 



AE e ~ 
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heavy neutrinos in the Galaxy are nonrelativistic, v ~ 300 km/ s is the velocity of neutrinos 
in the Galaxy) and from the continuum electron radiation resulting from decays of unstable 
particles. However, as we can see from Fig.^|(a) and Fig.|(a), due to energy losses the 
resulting electron spectrum becomes significantly softer with a sharp edge at E = m u . It 
appears that the electron spectrum cannot provide a clear signature for the heavy neutrino 
annihilation in the galactic halo. 



As was noted in |fL8| , 24 an important constraint on m u could be obtained by investigating 
the high energy part of the positron spectrum resulting from neutrino annihilation in the 
Galaxy. The flux of positrons from neutrino annihilation will be the same as in the case 
of electron production but the background conditions are significantly better. Up to now 
the positron flux has been measured only at energies < 30 GeV. To estimate the expected 
positron flux at higher energies we perfomed a 3-parameter fit to experimental data for the 
e + /(e + + e - ) ratio in the form 



R{E) = (el+jlS, d5) 

with three free parameters ai, ci2, a 3 . The result of our fit is presented in Fig.[|. In the 
energy region where separate measurements of electrons and positrons are available (see 
Fig.|l]) a significant excess of electrons has been found. At higher energies we can use either 
an extrapolation of our fit or the theoretical prediction. If we assume, for instance, the 
validity of the model of the dynamical halo, then the ratio e + /e~ is found to be less than 
0.02 at E e > 30 GeV [p3| . In this case we expect (see Fig.^|(b)) a clear signature for 
heavy neutrino annihilation by an edge in the positron spectrum, which corresponds to the 
annihilation process vv — > e + e~. The detection of an excess of high energy positrons and 
a sharp edge of the positron spectrum (at E = m v ) would be a clear signature of the 
annihilation of Dirac neutrinos in the galactic halo |18], , since the direct annihilation of 
massive Majorana neutrinos into light fermions in the Galaxy is severely by suppressed in 
the nonrelativistic case due to angular momentum conservation. 



For the case m v = 300 GeV the resulting electron spectrum is well below the background 
due to a sharp decrease of the neutrino number density when the neutrino mass increases 
above 100 GeV, n u ~ m~ 3 . 

The analysis of the p spectrum is similar to that of electrons and is based on expression @. 
However, it is important to note that in the case of antiprotons we can neglect their energy 
losses in the galactic disk due to mass suppression of synchrotron radiation and inverse 
Compton scattering in comparison with the case of electrons. Since antiprotons are trapped 
in the galactic magnetic field and have a rather large confinement time of the order 10 7 yr 
we can use the same expression (|9[) for an antiproton flux on the Earth which appears after 
a heavy neutrino annihilation and a sequential fragmentation of gluons and quarks. At low 
energies the ratio of antiprotons to protons in the cosmic rays is extremely low, namely less 
than about 10 -4 for E ~ 1 GeV (see Fig.|2|). Therefore, if a sufficient amount of antiprotons 
was created in the hadronization process the antiproton flux can excess the background and 
this could point at v v annihilation in the halo. 
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We present the results of our calculations of antiproton spectra in Figs.^(c), |(c). We ignored 
a modulation of the antiproton flux due to solar wind effects in the region above 1 GeV. We 



took the interstellar cosmic ray proton spectrum to be [25 



j; xp = 1.93 (^j E^ 7 cm~ 2 sr- 1 s- 1 GeV- 1 . (16) 
where, E^ in = E — M p is a kinetic energy of the proton. 

The dotted line shows our fit to the experimental data for the p/p ratio (see Fig.0) (in the 
form as Eq.flToD) multiplied by the proton flux ( |T6"D for the region above 13 GeV. As it is seen, 
the antiproton signal is too small to be detected in the case oivv annihilation. Even if we 
compare our results with standard leaky box model prediction for the pp ratio Rlb ~ 10~ 4 
in the region 10 — 100 GeV the results are not favourable and cannot impose any significant 
constraint on the model. 

Finally, we discuss the photon spectra resulting from v v annihilation. In the photon case 
there are no energy losses on the way to the Earth and the confinement time r in eq. ([|) is 
replaced by an integral over the line-of-sight which depends on the galactic latitude b and 
longitude / H, ||: 



jt = v -f^ n 'ji aI ^ ,l7) 



where for the diffuse 7 component factor I(b,l) ~ 1.24 (high galactic lattitudes) and a « 
10 kpc. The diffuse photon background is not known for E 1 > 1 GeV. Using data from the 



MeV region and extrapolating to higher energies, the photon spectrum is inferred to be [22 



J™ p = 8-10~ 7 £- 2 - 7 cm~ 2 sr' 1 s- 1 GeV' 1 . (18) 

The results of our 7 calculations are shown in FigfJ(d) and Fig.|](d). As we see in the case of 
m u = 100 GeV the neutrino annihilation leads to an indirect production of photons through 
particle decays above the expected background. However it seems that the continuous photon 
spectrum is not sufficient to give us a clear signature of the neutrino annihilation in the 
galactic halo. The difference between the predicted spectrum and the expected background 
is not very significant and the flux is small. 

We should note that our calculations contain significant uncertainties in astrophysical as- 
sumptions especially about the laws of propagation for p and e + in the Galaxy and the local 
halo density. 



5 Conclusion 



In this paper we have carried out a detailed analysis of fluxes of cosmic ray antiprotons, 
positrons, electrons and photons to be expected from the annihilation of relic heavy stable 
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neutrinos in the galactic halo. We used the standard electroweak model including one ad- 
ditional fermion family and the idea of neutrino condensation in the gravitational field of 
collapsing matter at the stage of Galaxy formation. Our conclusion is that it is difficult on 
the basis existing experimental data on cosmic ray spectra to constrain parameters of heavy 
neutrinos. However, it seems that future experimental study of cosmic rays and especially 
their e + and 7 components in the high energy region, could give important information on 
dark matter particles. In particular the detection of the anomalous positron output corre- 
sponding to the monochromatic annihilation line at the edge of the positron spectrum would 
be a clear indication for the annihilation of Dirac particles of dark matter in the galactic 
halo. Also the investigation of the continuum spectrum of diffuse 7 rays above 1 GeV could 
be used in order to constrain the parameters of dark matter particles or to provide a confir- 



mation of a signal seen in other experiments. As it was mention in [20] a possible strategy for 



future experiments is to combine positron detection with those for antiprotons and photons. 



Note that direct detection experiments p3| using semiconductor ionization detectors ex- 
cluded heavy Dirac neutrinos as a major component of our galactic halo in the mass region 
10 — 1000 GeV. However, the investigation of cosmic ray spectra could play an important 
role in the case of multicomponent dark matter and may be used for species (in particular, 
neutrinos) that give a negligible contribution to the total energy density of the Universe. 

We also emphasize the possibility of searching for very heavy neutrinos in accelerator ex- 



periments by the process e + e~ — > v v 7 [35[ which allows investigation in the region m v ~ 
rriz/2, where it is extremely difficult to get cosmological constraints on neutrino masses due 
to a very low neutrino number density in this mass region. If heavy neutrinos exist there 
could be also an interesting hadronless signature for Higgs meson strahlung production at 
accelerators e + e~ — > Z° H° —>■ Z° vv — > l + l~ vv and this mode could be the dominant one. 
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Table. 1. Total cross sections. 







= 100 GeV 




= 300 GeV 






= 3.90 GeV 




= 10.6 GeV 


Final state 


otou (nb) 


otou (nb) 


dd 




3.30 




0.245 


uu 




2.50 




0.189 


ss 




3.16 




0.245 


cc 




2.52 




0.192 


bb 




3.26 




0.287 


tt 








40.0 


e~e + 




0.737 




0.0554 


VeVe 




1.48 




0.110 






0.743 




0.0554 






1.48 




0.110 


T - T + 




0.724 




0.0571 


V T V T 




1.48 




0.110 


W + W~ 




16.2 




121 


£ over all channels 




36.2 




162 
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e + /(e + + e ) ratio Positron flux 

Figure 1: Measured data for the ratio e + /(e + + e~) (a) and for the positron flux (b) 
together with the curves resulting from the fit discussed in the text. 
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Figure 2: Measured data for the ratio p/p (a) and for the antiproton flux (b) together with 
the curves resulting from the fit discussed in the text. 
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Figure 3: Results of the simulation for m h 
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100 GeV . (a) — electrons; (b) — positrons; 



(c) — antiprotons and (d) — 7's. The histograms present the predicted particle fluxes for 
an energy resolution AE = 1 GeV. The dashed lines show the observable fluxes for present 
day experiments (see text). 
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Figure 4: Results of the simulation for m v = 300 GeV. (a) — electrons; (b) — positrons; 
(c) — antiprotons and (d) — 7's. The histograms present the predicted particle fluxes for 
an energy resolution AE = 1 GeV. The dashed lines show the observable fluxes for present 
day experiments (see text). 
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